Higher order corrections to the Newtonian potential 
in the Randall-Sundrum model 
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The general formalism for calculating the Newtonian potential in fine-tuned or critical Randall- 
Sundrum braneworlds is outlined. It is based on using the full tensor structure of the graviton 
propagator. This approach avoids the brane-bending effect arising from calculating the potential 
for a point source. For a single brane, this gives a clear understanding of the disputed overall factor 
4/3 entering the correction. The result can be written on a compact form which is evaluated to high 
accuracy for both short and large distances. 
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I. INTRODUCTION 

During the recent years, there has been a lot of ac- 
tivity on theories with extra dimensions. In contrast to 
the Kaluza-Klein picture where the radius of the extra 
dimension is usually taken to be the Planck length, these 
new theories allow for millimeter-sized extra dimensions 
as first proposed by Arkani-Hamed et al. or even 
an infinite extra dimension as proposed by Randall and 
Sundrum 2, 3J. These models were originally constructed 
to solve the hierarchy problem between the electroweak 
scale and the Planck scale and usually require Standard 
Model fields to be trapped on a 4-dimensional surface, 
called a 3-brane, to avoid conflict with experiments in 
particle physics where no signs of extra dimensions have 
been detected so far. Since gravity per definition can 
propagate through all dimensions, one important ques- 
tion is then whether gravity is still effectively localized 
on the brane, and Newton's law reproduced for distances 
larger than about a millimeter. 

Starting with Randall and Sundrum 0|, several other 
authors have calculated the Newtonian potential in the 
RS model, finding that V{r) - i[l-hA(r)], i.e. a 
dominant 4-dimcnsional term together with corrections. 
Although everyone agrees with the general form A ^ 
l//i^r^ of the leading order correction (for large dis- 
tances), there seems to be an ambiguity concerning the 
numerical coefficient. By looking at the scalar part of 
the graviton propagator, the result A = l/(2^^r^) is ob- 
tained 0, Q , whereas calculating the component /iqo of 
the metric perturbation due to a matter source on the 
brane gives the result A = 2/ {i^JL'^r'^) US 110. In the 
latter calculation, complications arise because the posi- 
tion of the brane can no longer be fixed in the fifth di- 
mension when introducing a matter source. The relative 
factor 4/3 between the above results is then considered 
to be the net result of this brane-bending effect. 
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In this paper we calculate the Newtonian potential 
from the graviton propagator for a critical brane (with 
cosmological constant A = and curvature k — Q), and 
show that the factor 4/3 is obtained by using the full ten- 
sor structure of the propagator instead of just the scalar 
part, without the need to include any brane-bending ef- 
fects since we don't have a matter source 0| . This also 
gives the same general expression for the correction as 
obtained by Chung et al. but derived in more detail 
here. Throughout this paper, we try to include most of 
the details in each calculation, hopefully making it pos- 
sible to follow without much prior knowledge on the sub- 
ject. The only parts left out are either straightforward 
or easily found elsewhere in the literature. 

Other works which are closely related to this paper 
include Q and ^ , where an induced 4-dimensional cur- 
vature term (i.e. a term ~ J d'^x^JgR in the action) 
modifies the Newtonian potential so that it behaves 4- 
dimensional even at very short distances, and [l^ . where 
n compact extra dimensions are included in addition to 
the large one. In the conditions for localization of 
gravity are addressed more systematically, and in it 
is shown that gravity is localized also for a non-critical 
brane, meaning that localization is a local property of 
the model, unaffected by the geometry of space far from 
the brane. The localization of other spin fields (spin 0, 
1, i and |) is studied in 15]. 

This paper is organized as follows: In section^we be- 
gin by solving Einstein's equations for the 5-dimensional 
RS model, thus obtaining its metric. In section Ulll we 
make a perturbation to the metric and derive the wave 
equation the perturbation has to satisfy, starting with a 
completely general metric and then specializing to the 
RS case. The wave equation is then solved, giving us 
the graviton wavefunction. From the wavefunction we 
obtain the graviton propagator in section Hvl and finally 
the Newtonian potential between two point masses on 
the brane, expressed as an integral over the mass ni of 
the gravitons. In order to obtain the correct normaliza- 
tion for the wavefunction, and also the correct integration 
measure over the mass m, we use a regulator brane at a 
distance y,- from the physical brane, and take the limit 
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Ur ^ oo at the end of the calculation. We also expand the 
potential in powers of the distance r between the point 
masses, obtaining the leading order corrections 2/{'in^r'^) 
and 4/(37r/xr) for large and short distances, respectively, 
along with higher order corrections. 



II. THE RANDALL-SUNDRUM METRIC 

We start out by assuming that the 5-dimensional met- 
ric has the form 
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ds^ = n^it, y)df - a^{t, y)j,jdx'dx^ - dy^, (1) and n{t, y) ^ a/do. H ere M is the 5-dimensional Planck 

mass, fi — -y/IAsl/G, and oq = a(t, y = 0) denotes the 
.^jjgj.g scale factor on the brane. The latter satisfies an equation 

similar to Friedmann's first equation 
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is the spatial metric. The metric is thus parametrized 
by two scale factors n(f, y) and a(<, y), with y being the 
coordinate in the fifth dimension, and y = taken as the 
position of the brane. We use the standard convention 
n{t,y = 0) = 1, which means that t is the proper time 
on the brane. The energy momentum tensor is taken to 
be Tmn 



rpB 



T, 



MN 



^MN, where 
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where J7 is a constant of integration, often called the dark 
radiation term for obvious reasons. In addition, we have 
the energy conservation equation 



Pb 



3— (P6 



Pb), 



(8) 



just as in standard cosmology. Note, in particular, the 
absolute value \y\ in ©, which is the result of the bound- 
ary conditions 



is the contribution from the bulk space (index B), and 

Tmn = diag(pb"-^Pba^7^J,0)(5(?/) (4) 

is the contribution from the brane (index 6). Finding 
the Einstein tensor and solving Einstein's equation is 
straightforward and has been done by many authors (see 
e.g. 16, 17, 18, 19]), so we simply state the results here. 
First, the Einstein tensor is given by 
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where a dot denotes the derivative with respect to i, 
and a prime the derivative with respect to y. Assuming 
that the bulk space only contains a cosmological con- 
stant, i.e. Pb = M^Ab = —PB, Einstein's equation 



1^ = -lAr 



Pb, 



^ - iM-3(2pb + 3pb), (9) 
no 



imposed by the delta function in T^j^. Here [a'] means 
the jump discontinuity in a' across the brane, i.e. [a'] = 
a'(y = 0+) - a'{y = 0") = 2a'(0+). If we write the 
energy density p& on the brane as pb — X + p where A is 
a constant, iQ takes the more recognizable form 



P 



3M2i 



A 
"3 



U 



(10) 



where we have identified the effective 4-dimensional 
Planck mass 



Mpi = {M-^x/ey^^^ = M^ye/A, 

and the effective cosmological constant 



i(Mp7A + AB) 



(11) 



(12) 



on the brane. The constant A is often called the brane 
tension, and is assumed to be positive. Clearly, re- 
produces standard 4-dimensional cosmology when U is 
"small" and A "large". We will not discuss more pre- 
cisely what this means other than saying that A has to 
be (much) larger than the typical densities /0„uci. of mat- 
ter and radiation during nucleosynthesis. 
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For the remainder of this paper, we will make the sim- 
plifying assumptions that pf, = A (the brane contains 
no components other than its tension) and that U = 0. 
Under these assumptions, the scale factor a{t,y) can be 
factorized as a{t, y) = ao{t)A{y), where A{y) follows from 



metric 



GAP' 



My) = { 



cosh /iy 



A 



■sinh^lyl. As < 0, (13) 



cos^y - ^^^3 singly], Aij > 0. 



It also follows that n{t,y) = A{y), and the metric is 
simplified to 



ds^ = A^{y)gf_,i,{x)dx^'dx'' - dy^, 



(14) 



where g^i,j{x) — dt^ — a'^{t)^ijdx'^dx^ is the stan- 
dard Friedmann-Robertson- Walker (FRW) metric on the 
brane. 

Our main focus in this paper will be the case of a 
vanishing effective cosmological constant, A = 0, which 
is often called a critical brane. We then get Kb = 
— iAf~^A^ < 0, meaning that the bulk space is AdS^, 

and A{y) — e^^'^L If we also assume a spatially flat 
universe, fc = 0, the metric becomes 



(15) 



This is the case originally studied in 0, 0] . One should 
note that A^ < is in fact the only possibility compatible 
with observations. Since the observed value of A is very 
small, the requirement that A S> Pnuci. obviously gives 
Mpi^A > |A|. From (O we therefore get Ab = 2A - 
Mpi^A « -Mpi^A < 0. 



III. LOCALIZATION OF GRAVITY 

The main goal of this paper is to calculate the correc- 
tion to Newton's law of gravitation due to the presence 
of the fifth dimension. Since the Newtonian potential 
is essentially the low energy limit of the graviton prop- 
agator, we must find the graviton wavefunction and in 
particular study its localization on the brane. The gravi- 
ton is described by the traceless transverse component 
hij of the spatial metric perturbation. The question of 
finding the linearized wave equation satisfied by hij has 
been addressed by many authors ^,20^ 21], but since 
this is an integral part of the calculation, we include the 
details here. 



A. General linearized gravity 

To begin with, we assume nothing about the back- 
ground metric g^j^, and consider the general perturbed 



(16) 



where also the number of dimensions is arbitrary. In the 
following, we will only keep terms to the first order in h 
and its derivatives. The inverse metric is therefore 



and the Christoffel symbols 
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From this, we see that the perturbation S^^^ to the 
Christoffel symbols is in fact a tensor. (The symbol 
; means covariant derivative with respect to the back- 
ground metric g^ ^.) Next, we find the Riemann tensor 

i/a/3 ~ i/a/3 ~^ i//3,a '-' va.p 

lyP CM I rt^ qP —VP — r'^ ^P 

= ^ta/3 + ^^ul3;a ~ ^ta;f3 ' (l^) 

and thus the Ricci tensor 



p pa p I na not 

J^Hf — f_i_au ^ ^i^v fiv;a pa;u 



, (20) 



where the last expression is obtained after a short and 
straightforward calculation. Finally, we obtain the cur- 
vature scalar 

R = ^R- h'^'R^u + h^'-^.. ~ ah , (21) 

and the Einstein tensor 

— -|- 2 ihoL^\y H" haiy;fi ^fj.v^a 

-^Rh^, + (h'^^Rap + ah- h''^,^fi) . (22) 

This can be somewhat simplified by using the trace- 
reversed components h^u = h^^, — ^g^^h, yielding 



2(^0/1;!^ ~^ ^au-fi ^ a^iiu Rhfii, 



(23) 



Throughout this paper □ always means the scalar 
d'Alembertian, i.e. □ = {1/ ^/g)^f^{y/gg'^''^I,), whereas 
= g^^^a^fs means the covariant Laplacian which 
acts differently on tensors of different rank. It is there- 
fore only for a scalar field (j) that V^c^ = Ocj). 

The invarianse of (|23|l under an arbitrary coordinate 
transformation can be used to simplify the expression. 
We choose to work in the gauge where h = and h^'^.j^ = 
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0. In order to utilize the latter condition in (|23|l . we use 

the general tensor identity 

Tal3;fj,i/ ~ Tafj-^y^ = TafjR p^^^ + T^^pR , (24) 

which can be proven by " brute force" by writing out the 
covariant derivatives explicitely. From this it follows that 

^Q/j;i/ ~ ^ fj.-.aiy + ^ fii^a + ^ ^Rav i (25) 

and H23|l is then reduced to 

-Vlh^, - Rh^, + 9^,,h''^Rc.p) . (26) 

Taking the trace, the first three terms are seen to cancel, 
yielding 

5E^'^, = li-Dh -Rh + Dh^f^Ra^p) = y/i"^i?a/3, (27) 

where D is the number of spacetime dimensions, and we 
have used the traceless condition h ^ 0. From now on, 
we can also ignore the bar, since h^^, = 

Einstein's equation for the perturbed metric can 
be written E^^ = M'^^ T^^, where Md is the D- 
dimensional Planck mass, and T^^ the perturbed energy 
momentum tensor. Assuming that both energy den- 
sity, pressure and " Z?- velocity" u^^ is unperturbed, i.e. 
5p ~ 5p ~ Q and 5u^^ = 0, we get 

Tt_,v = {p + p)u^u„ - pg^i, 

^ {p + p)u^u„-p{g^^ + hf_,^), (28) 

meaning that ST^^^, — —ph^j^^. The background metric 
is determined from the unperturbed energy momentum 
tensor, i.e. E^^, — T^^^, and the equation for the 

perturbed Einstein tensor is therefore 

5E^, = Ml-"" 5T^, = -Ml-'^ph^, . (29) 

From (|27|l we then get 

5E% = ^h^^R^p = -M^-'^ph'^^ = 0, (30) 

since h is traceless. When inserted into (|26(1 . this gives 

-S/lh^, - Rhf,, + 2Ml-^ph^,, = . (31) 



The first four terms are often called the de Rham- 
Lichnerowicz operator A applied to /i^^: 



Ahf,, = 2/i",i?"^,„ + h%,Ra. + h\R,^ - Vlh^, . (32) 



The general equation (|31|l can be applied to any geometry 
we are interested in for tensor fluctuations. 



B. Linearized gravity in the RS model 

We will now apply 13111 to the RS background metric 
JH). In addition to h^'^'^ -n = and h^'\i = 0, we choose 
the gauge where /i^4 = ^144 = 0. Having only one brane, 
this was shown possible in a more detailed discussion 
about gauge transformations by Boos et al. The 
condition that ^144 = corresponds to having no massless 
scalar field (radion) in the problem. (With two branes 
such a gauge choice is not possible, and the radion will 
contribute in that case.) 

We now have a total of 10 gauge conditions, reducing 
the number of independent components of the symmet- 
ric tensor Hmn from 15 to 5. These components must 
then be separated into scalar, vector and tensor parts 
with respect to the spatial metric 7^ . Here we choose 
to work only with the tensor part which describes pure 
gravitational waves, meaning that we set h^o = 0. The 
remaining part hij therefore have two independent de- 
grees of freedom. If we also factorize out by letting 
hij o?hij, the full perturbed metric can be written 



ds^ = n^{t, y)dt^ — a^{t, y)ijij + hij)dx^dx^ — dy^ . (33) 



The indices of hij are now raised and lowered by 7^, 
and the gauge conditions h^^ -n = and h'^ m = 
are reduced to '^jh^-^ = and /i\ = j'^^hij = 0, where 
Vj is the spatial covariant derivative. Thus hij may be 
regarded effectively as a 3-dimensional tensor, with t and 
y as additional parameters. The different terms in l|31|l 
are then easily found using Mathematica |2^], with the 
result 
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(34) 



The last equation also follows directly from the Einstein 
tensor since Eij = M'^T^j = M'^pa^^ij, with p = 
Pb + PbS{y). Adding everything together, we see that 
most of the terms cancel, and we are left with 



From we see that the x-dependence of /ly can be 
separated from the t- and y-dependence. Writing hij = 
i/;(i,2/)Gy(x), (123 is solved by 7'='VfeV,Gy = -a^Gy 
and 
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A tedious but straightforward calculation shows that the 
5-dimensional Laplacian can be expanded as 
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do 
2fc4 



ij{t,y) = 0. (39) 



The eigenvalues a' of the 3-dimensional Laplacian can be 
found in but we will not need them here. Instead, 
since we are considering the factorized metric H14|l . it is 
more convenient to separate the y-dependence only. 

Substituting a{t,y) = ao{t)A{y) and n{t,y) = A{y) 
into (|37|l . it simplifies to 



2k 
a, 



with the solution hij(x,y) — Gij{x)^{y), where 



h,j = , (40) 



(Removing the y-dependence and setting n = 1, we re- 
cover the result for the 4-dimensional FRW metric as 
obtained by Ford and Parker HJ].) We should compare 
(|^ to the wave equation of a scalar field with mass 
in five dimensions: 



(□ + m2)0 = 
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do 



2 7'''VfcV/ 



= 0. 



(38) 



The only difference between the two expressions lies in 
the tensor structure of the term 7'^'VfeV;/iij. The gravi- 
ton field has also acquired an effective mass nil = 2k /a'^, 
which at the present epoch is of the order 10"'^^ eV 
(at most), and therefore totally negligible at short (i.e. 
non-cosmological) distances. For a spatially flat universe 
{k — 0), l|T7)l is reduced to Dhij = 0. The components 
of hij independently satisfy the Klein-Gordon equation 
in this case, and the graviton field is thus equivalent to 
a set of two independent scalar fields. 
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52 -t- -^do - —I'^Vkyi + — +m'\ Gy (x) = 0, (41) 
ao a^ a^ J 

4- A' 777^ 

<^"{y) + —^'{y) + ^m = 0, (42) 



and m is the eigenvalue following from the separation of 
variables. This is the same result as obtained by Brevik 
et al. 01 . Explicit solutions to (|^ can be found by 
separating the t-dependence, but we will not need them 
in the following. The mass spectrum will be determined 
from (02). 

One should also note that the result (|42|) can be ob- 
tained without demanding that ft.^o ~ when the 5- 
dimensional metric has the simple form (|14|l . Instead, 
the perturbed metric can then be written 



dy' 



(43) 



with the gauge conditions h'^" -y — and h = g'^'^h^i, — 0, 
where the indices of h^^ are raised and lowered by g^^. 
Thus h^i, is effectively a 4-dimensional tensor, with five 
independent degrees of freedom. Going through the same 
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calculations as above, starting from H31|) . we then get the 
result 
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l2 



2A 



h^u^O, (44) 



where is the 4-dimensional Laplacian. The appar- 
ent difference in the effective mass term compared to 
H4()|l is simply because now we have not factorized out 
ao(t) from the perturbation. Writing h^^, — a'^{t)h^i, 
and using Friedmann's equation l(Tn|l . h^^ can be seen 
to satisfy (|^ . provided we also set h^Q = in or- 
der to expand V^. Equation 144|l is solved by setting 
hfj,u{x,y) = Gf^^{x)^{y), where $(2/) satisfies (|^ . and 



2A 
T 



0. 



(45) 



Thus, TO can be regarded as the effective mass of the 
5-dimensional graviton field as observed in four dimen- 
sions, with a small correction from the cosmological con- 
stant. When A > 0, even the m = state is massive. 
This appears to create a serious problem, since massive 
gravitons (in 4 dimensions) have five polarization states 
whereas massless gravitons only have two polarization 
states, and all physical results should be continuous in 
the limit A ^ 0. This problem and its solution has been 
discussed by Karch and Randall . 



C. The graviton wavefunction 

Next, we proceed to solve (j^ . Since this has already 
been done we will only sketch the derivation. 

First, we simplify the equation by making the substitu- 
tion y — > z{y), where 



(46) 



and by writing $(y) as 

$(y) ^ A-'^''^u{z). (47) 
This leads to the Schrodinger like equation for u{z) 



[-dl + V{z)\ u{z)^rr?u{z), 
where the effective potential V{z) is given by 

V{z)^\{A'f + lAA". 

Notice that the line element (|14f) takes the form 
ds^ = A^{z) (g^.^dxt^dx" - dz^) 



(48) 



(50) 



when expressed using the coordinate z. Therefore z is 
often called the conformal coordinate. 



Focusing on the critical case A 
g-Mlyl^ we get the result 
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where A(y) 



z sgn(y)- I 
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(51) 



where the delta function originates from the double 
derivative of \y\, and we have used the general relation 
S{z - z') = A{y')6{y - y') together with A{0) = 1. The 
potential for the case A > can be found in e.g. 18, 2^, 
and in |Q| the case A < is also included. The delta 
function in the result means that u{z) must depend on 
the absolute value \z\ (just as the delta function in Tmn 
resulted in a{t,y) depending on \y\). Matching the delta 
function in V{z) with the delta function from u"{z), we 
get the boundary condition 



2u'(0) + 3//u(0) = 0, 



and the equation for z > 
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u{z) = 0. 



(52) 



(53) 



The general solution to (|53|l is easily shown to be u{z) — 

^^rTT^^^AJ2[fil + flz)]+BY2[f{^ + ^lz)]'j for arbi- 

trary constants A and B, where J2{x) and Y2{x) are the 
Bessel functions of order two of the first and second kind, 
respectively. In order to satisfy the boundary condition 
(|52|l . we must choose the linear combination 



xlY2[^{\ + M 



-^J2[^(l + Mkl)][, (54) 



where N^a is a normalization constant (to be found later). 
This expression is valid for all m, so we have a continuous 
mass spectrum for m > 0. The solution for m — Q can 
be found either directly from H53|l or by taking the limit 
TO ^ in (|54|1 . and is given by 



u^{z) = NQ{l + ^l\z\)-^'^. 



(55) 



Thus gravity is localized on the brane, since the massless, 
bound state U{){z) can be normalized, and has a sharp 
peak at z = 0. As we will see later, the massive states 
(|54|l give the deviation from Newton's law. Solutions for 
the case A > are given in ^ J8„ ,2^ , and for A < in 
(49) [l|. 



IV. CORRECTIONS TO NEWTON'S LAW 

Having found the graviton wavefunction, we will use 
this result to calculate the static gravitational potential 
between two point masses both located on the brane. As 
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is well known from field theory, the (Fourier transformed) 
potential is obtained by looking at a scattering diagram 
where two particles interact through the exchange of a 
virtual graviton, in the limit where the energy of the 
graviton goes to zero: 



F(k) = hm 

fct>-»0 




(56) 



We therefore need to know the 5-dimensional graviton 
propagator, as well as the interaction between a gravi- 
ton and a massive particle on the brane. The propagator 
must in principle be found by inverting the quadratic 
terms in the effective graviton Lagrangian [22l |. How- 
ever, since we have expressed the 5-dimensional gravi- 
ton wavefunction directly in terms of the 4-dimensional 
ones through hf,^{x,y) = J2m^]^iy(^)^('^^y)^ know 
that the 5-dimensional propagator must be similarily ex- 
pressed by the 4-dimensional propagator, and can there- 
fore immediately write down the answer: 
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(5) 



p{x,y;x\y') = {0\Th^^{x,y)haf3{x' ,y')\0) 



E 



$(m,y)<i>*(m,y')<r/](^-^'), 

(57) 



where ^^t''™/3(^' propagator of a 4-dimensional 

spin-2 particle with mass m. Taking y ^ y' — and 
using $(m, y) = A{y)^^^^u{m, z), we get 



Ko)i^<j:j(x,x') 



(58) 



Everything is now happening in just four dimensions. 
The only trace left of the motion of the graviton 
through the fifth dimension is that we have a tower 
of massive 4-dimensional gravitons with a non-trivial 
"weight" |u(m, 0)p. More importantly, the propagator 

^livafii^' a;') is completely described in the 4-dimensional 
FRW space with metric g^jy, where standard results can 
be applied directly without worrying about the fifth di- 
mension. 

Focusing again on the critical case k = and A = 0, 
the 4-dimensional space is flat (g^i^ = with the 

result 



d*k 



(27r)4 fc2 _ ^2 + ie 



-ik-{x—x ) 



(59) 



where the polarization tensor can be chosen as (see e.g. 
[25| for more details) 



1 



{•ntiaripf} + rj^p-qva - ■q^vriaf}) , (60) 



1 



2m? 
1 



kp + ri^ipkyka + rjjjak^kp + rj^pk^^ka) 



6 V 



k^ ky 



Vap H J 



kakp 



(61) 



Next, we consider the interaction. For a Lagrangian 
C in D dimensions the action is S* = / d^Xy/gC. The 
energy momentum tensor Tmn corresponding to the La- 
grangian C is defined by varying the action with respect 
to the metric: 



SS : 



J d^Xy/g^TMN 



Sg 



MN 



(62) 



Since Sg^'^^ = — /i*^^ to the lowest order in /i, we there- 
fore get the interaction Lagrangian 



iMNiT' 



(63) 



So far, we have considered the metric perturbation Hmn 
to be dimensionless, since gMN is dimensionless. But 
when treating it as a field of particles (i.e. gravitons), this 
can no longer be the case if we want to have a canonical 
kinetic term in the Lagrangian. Instead, Hmn has to have 
a dimension dim[ft,] = [D — 2)/2 in D spacetime dimen- 
sions. Since /imjv is a 5-dimensional field, we must there- 
fore have dim[/i] = 3/2. From h^^[x,y) = G^v{x)'^{y), 
we then get dim[$] = dim[u] = 1/2, since G^i, is a 4- 
dimensional field and should therefore have dimension 
1. This is also required in order to make equations 1571 - 
I61|l consistent. Again, since h^^ is 5-dimcnsional, we 
should use the 5-dimcnsional Planck mass M and let 
ht^'^ ^ in (ggi. Finally, using = /iM^j 

which follows from and H12|) . we get the result 




jjLV = 



y/JIMpi 



(64) 



Before proceeding, lets give a little thought to the as- 
sumption A = and fc = 0. For nonzero A, we should ex- 
pect corrections to Newton's law of the form (lAjr^)" for 
some positive constant n. Because of the very small ob- 
served value of A, this correction will only be noticable at 
very large distances. But the Newtonian potential is only 
a meaningful quantity at short (i.e. non-cosmological) 
distances, so for our purpose the approximation A « 
should be perfectly valid. The same is also true for a 
potentially non- vanishing fc. From a pure mathematical 
point of view, however, it would be interesting to con- 
sider an arbitrary value of A. But in that case, one has 
to use the full deSitter propagator instead of 

Another way of looking at this problem, is to note that 
the effect of the fifth dimension and the effect of A and 
k takes place at completely different scales. The effect of 
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the fifth dimension dominates at very short distances (as 
we will see), where we can set A = and fc = 0, whereas 
the effect of A and k dominates at very large distances, 
where we can ignore the fifth dimension altogether and 
simply use 4-dimensional FRW space. Trying to include 
both effects in the same calculation will only make things 
unnecessarily complicated. 



However, the second integral as it stands is divergent for 
all m and m', not just for m = m', and therefore has to 
be regulated in some way. As mentioned in this can 
be done by introducing a regulator brane at some large 
but finite position Zr, and then taking the limit — > oo 
at the end. Since this point seems to be taken very lightly 
in the literature, we will go through the calculations in 
detail here. 



A. The Newtonian potential 

Combining the above results H59|l and (|64(l for the prop- 
agator and the vertex and summing over m, we finally 
obtain the gravitational potential 



y(k)= ^|u(m,0)p 



t ^1 ^f_iuafS-^2 



(rn) rpCtP 



IP 



(65) 



fcO— 



Taking both particles to be point particles at rest with 
masses mi and m2, the energy momentum tensors can be 
written rf''(x) = mi(5(x)M''u'^ =?> T|"'(k) = miw'^w" = 
TOiJqJq, and likewise for T2^. We therefore see that 
only the 0000-componcnt of the polarization tensor con- 
tributes to the potential: 



V{k) 



8TrGmim2 



J2\u{m,0)\ 



2 ^00 00 V 
k2 



(66) 



(Here we have written Afpj^ = 8nG, where G is the grav- 
itational constant.) Now, since fco — 0, the polarization 
tensor is easily found from 16U|I and (|61|l : 



d('") 
0000 



(k) = 




m — , 
m > , 



(67) 



with the result 

87rGmim2 



V^(k) 



lju(0,0)| = 



|M(m,0)|^ 



k2 



3 ^ k2 

m>0 



(68) 



or 



Gmim2 

jJLT 



l"(0,0)p + l ^ |u(m,0)pe— 

m>0 



(69) 

Thus, we see the importance of using the full tensor struc- 
ture of the propagator, instead of just the scalar part 
which is commonly used. Without the tensor part we 
would miss the relative factor 4/3 between the massless 
and the massive modes. 

In order to proceed further, we need the normalization 
constants Nm and Nq in (|54|l and H55|l . which in principle 
are found by requiring 



|u(0,z)pdz = 1, 
u(m, z)*u{m' , z)dz = S{m — m') . (70) 



B. Using a regulator brane 

With a second brane located at y = yr, we still assume 
that the points y and — y are equal, which means that 
the fifth dimension is now compactified on an orbifold 
S^/Z2. The two branes thus represent the endpoints of 
the fifth dimension. Einstein's equation in the bulk space 
and the boundary conditions on the physical brane are 
unchanged by the presence of the regulator brane, and 
the only way it manifests itself is as an additional set of 
boundary conditions. (Several works have been done in 
the two brane scenario where both branes are considered 
physical, see e.g. . In this paper we are only using the 
second brane as a regulator, therefore ignoring its cosmo- 
logical implications.) From the energy momentum tensor 
'^MN = diag(/9rn.2,pra^7y , 0)(5(?/ - yr) on the regulator 
brane, we get 



a 
n 



v=yr 



v=yr 



iM-3(2p, 



3pr), 



(71) 



where [a']y=y^ = a'{y = y+) - a'{y = y^) ^ -"^a'iy^ ). 
Since a{y) is fixed by I© or H13fl . the energy density pr 
on the regulator brane can not be chosen independently, 
but is in fact directly related to the energy density pb on 
the physical brane. 

The condition (|71|l gives an additional term in the po- 
tential for the wavefunction, 



V{z)^ 



3/i(5(z) 



3^1 



6iz-zr), (72) 



4(l-h/i|z|)2 ' 1+pZr 

and thus the new boundary condition 

2u'{Zr) — 3A' {yr)u{Zr) — 2u' {Zr) + — u{Zr) — , 

1 + flZr 

(73) 

in addition to (|52|l . The zero mode H55|l is trivially seen to 
satisfy this new condition. But the continuous spectrum 
of massive modes is reduced to a discreet spectrum. Us- 
ing H54|) together with well known identities for the Bessel 
functions, H73|) can be simplified to 



Mfil+pZr)] 



ia.n{mzr - ^) , (74) 
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which means that the mass is approximately quantized 
in units of it / Zj.: 



constant cancel, as it should. We then finally obtain the 
gravitational potential 



rrir, 



■, 71 = 0,1,2,3,.. 



(75) 



The transition from a sum over to to a continuous integral 
in the limit Zr oo is therefore given by 

Z /"oo ^ 

/ /("i) — / /(to)— Am-^ / f{m)—dm. (76) 



Next, we consider the normalization constant Nr, 
when using a regulator brane, which is found from 



Ui{z)uj{z)dz = 2 / Ui{z)uj{z)dz — 6ij , (77) 



since the eigenfunctions Ui{z) = u{mi, z) in 154|) are real. 
When i ^ can use (I48II and write the integrand as 



Ui{z)uj{z) 



7 7 

m- — m- 



h(z)<(z)-u,(z)<(z)]. (78) 



Performing a partial integration and using the boundary 
conditions at both cndpoints, the integral then vanishes 
identically in the case i ^ j- For i — j we get 




u {m,z)dz— / N^{1 + fiz) 
Jo 

x{Y,[fi^ + ^,z)]^^^^Mfi^ + ^iz)]^ dz. (79) 

In the limit z,. — )■ oo we can use the asymptotic expres- 
sions for I2 and J2 for large arguments, with the result 

(..iv^)-i- hm 1 

Zr^oo Zr Jo """^ 

f Yii-) V 

X |sin(TOz - ^) - J^JW^ cos(toz - x) | '^^ 



2^ 

TTTO 



2lJ.Zr 



For the zero mode (|55|l we take the limit Zr 
with the result 



^0 



dz 



1 -1 



-00 (1 + Mkl)^ 



(80) 

(81) 
cx) directly, 

(82) 



Inserting everything into Ht)9|l , we see that the factor z^ 
from the integration measure and from the normalization 



Vir) 



GTOITOT ( r, , , 4 ry , ^ ^^Zr , 

—{u^(0,0) + - u2(TO,0)e"™'^— dm 

nr [ 3 Jo TT 

GTO1TO2 



Ji(^)r2(^)-ri(^)J2(^) 



-dm 



GTO1TO2 



e dm 



(83) 



where we have used the general identity Jn{x)Yn+i{x) — 
Yn{x)Jn+i{x) = —2/Trx in the last line. The first term 
clearly gives 4-dimensional gravity, and the integral over 
the massive modes gives the correction due to the fifth 
dimension. Since the only parameter in the integral is /ir 
(after substituting the dimensionless variable n — m/^), 
we see that the scale where the correction becomes im- 
portant is given by the curvature radius /i^^ of the fifth 
dimension. For distances r much larger than this, the 
exponential factor ensures that 4-dimensional gravity is 
restored. The result H83|l is the same as that obtained by 
Chung et al. Q. 

The integral in (|83|l must in general be done numer- 
ically, which is straightforward using Mathematica, and 
the result is shown in figure ^ It is also possible to 
expand the result in powers of /xr for the two regions 
/ir ^ 1 and iir 1. Writing the full potential as 
V{r) — Vo(r')(l + A), the relative correction A has the 
series expansion 



- /ir In /ir 



4 1 1 
37r/ir 3 27:' 

-H0.089237810/ir+C'(/iV2), /ir < 1, 
2 41n/ir 16 -12 In 2 (84) 



3//^r^ 



(In jirY 



3/i4r4 



/ir > 1. 



A detailed derivation of this result is included in ap- 
pendix ^ Note in particular the leading order term 
A ~ 4/(37r/ir) for short distances. When /ir <C 1, 
A ^ 1 and the entire gravitational potential behaves 
like V ~ 1/r^. Gravity is therefore 5-dimensional at short 
distances. This should not come as a surprise, since when 
the distance r is small compared to the curvature radius 
of the fifth dimension, spacetime looks almost fiat 
and we should get the same result as with 5-dimensional 
Minkowski space. 
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log A 
4 

2 - 




log fir 



FIG. 1; The figure shows the relative correction A in 1)83^ . 
where V{r) = Vo(t')(1 + A), as obtained by a numerical 
integration using Mathematica. The leading order terms 
A ~ 4/(37r/xr) (short distances) and A ~ 2/{Sfi^r^) (large 
distances) are clearly seen from the figure, and the transition 
between the two regions occurs around /ir ~ 1. 



V. SUMMARY 

In this paper, we have gone through the derivation of 
the Newtonian gravitational potential in the RS model, 
from start to finish. We have concentrated on the critical 
case (A = and fc = 0) , and have obtained the result (for 
large distances) V ^ + 371^ + • . •) from the graviton 
propagator, by including its full tensor structure. This 
is exactly the same result as what is found from a mat- 
ter source on the brane when the brane-bending effect 
is included. As for the non-critical case, since both A 
and /c/flg are very small (from observational constraints), 
corrections due to A and k will be noticable only at cos- 
mological distances, and the correction due to the fifth 
dimension must therefore be the same as for the critical 
case 

Throughout this paper, we have assumed that the dark 
radiation term vanishes (U = 0), and we have ignored 
contributions from matter and radiation on the brane by 
setting p = p„i + pr = 0. In a realistic model, how- 
ever, these terms have to be included. Because of the de- 
generacy between dark radiation and ordinary radiation, 
the two of them can not be distinguished by studying 
cosmological evolution alone. But since no observation 
has been made that requires the presence of dark radia- 
tion, we can probable assume that it is small, i.e. that 
\U\ < Prat/M^,. 

It is interesting to note, however, that the scale factor 
a(i, y) may change dramatically far from the brane when 
including even a tiny value of p or U , which is seen from 
the following example: For As < and A = we assume 
that all parameters of the model (M, p, and A) are of the 
order of the Planck scale, which means that we ignore 
the p^ term in H10|l . Combining this with ijHJl, the scale 
factor reduces to 



g'ft, y) 



-2Miyr 



PI 



U sinh py 

p'^al 



(85) 



Thus, both p and U introduce a horizon at a finite dis- 
tance from the brane, unless U is negative, in which case 
the scale factor reaches a minimum and then diverges as 
y ^ 00. The same two effects arise when A > and 
A < 0, resp ectively. However, as shown by Karch and 
Randall gravity can still be localized in these cases, 
since localization only depends on the scale factor close 
to the brane. We should therefore expect that matter 
and (dark) radiation only play a significant role in the 
very early universe, when p and U/gq are of the order of 
the Planck scale. 
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APPENDIX A: SERIES EXPANSION OF A 

We are considering the integral 

(Al) 

where we are using the dimensionless integration variable 
n — m/p. In the following, we will need the series ex- 
pansion of f{n) for small and large values of n, which is 
found using Mathematica: 



37r2 7o J2(n)+Yj2(„) n 



2 2 
/(n< 1) =3^+3 



1 



7 — In 2 + In rt ) n 



-0.107720074417+ -(Inn)^ 



7 — In 2 In n 

' 12 ' 



fin » 1) 



1 



21 



633 



37r 27rn2 32?™* 256Tm^ 



+ 



(A2) 
(A3) 



1. Short distances, pr <^ 1 

In the limit of short distances, the exponential factor 
g-fj-rn hecomes important when n is large and we can use 
(|A3p . However, this series can not be integrated term 
by term from to cx3, since all terms (except the first 
one) diverges in the limit n ^ 0. We therefore split 
the integral into two parts, ~ Iq ~^ //v°' '^t^^''^ N is 
assumed to be large. In the first part we expand e^^''", 
and the second part is integrated term by term, using the 
general formula 



00 — /^rn 



fe-2 



N n- 



1=0 



(fc-1)! iV^' 



l-i 



(A4) 
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Here T{a, z) is the incomplete gamma function r(a, z) = 
IT t°'~^G~^dt. Writing the expansion of f{ri) as f{n) — 
Sfe^o C'2fc"~^'^, we then obtain 



^ ' I n''f{n)dn ' 



k=0 



kl 



37r/ir 



E 

2fc-l 



(2fc-l)! Ar2fc-i-» 



(2fc- 1)! 



T{0,firN) 



(A5) 



and we can then extract and calculate the coefficient in 
front of each power of fir in turn. From HA6|I . we see that 
we also get several logarithmic terms. For each power of 
fir, we take the limit TV — > cx), therefore ignoring all terms 
with a negative power of N. Using the notation 



A= ^ Afe(/zr)'^ + logarithmic terms, (A7) 



This expression can now be expanded in powers of fir, 
using 



r(0, firN) = -7 - In - In TV - ^ 



k=l 



(-firNf 
kl ■ k 



we first find rather trivially that A^i = 4/ (Stt), and then, 



, (A6) using numerical integration: 



A, 



lim 

N^oo 



lim 



/ f{n)dn- ^ = / 
Jo 7o 









fin) ^ 


Jo 


6n 



dn = — , 
3 



N 



2N^ 1 
nf(n)dn H IniV 



1-7 
27r 



fin) 



3tt 



dn 



fin) 



1 



Btt 27rn2 



dn 



1-7 
27r 



0.08923780957038536 , 



Ao = 



lim 



n j(n)dn 



lim 

N^oo 



1 

3! Jo 



N 



n^ f(n)dn 



97r 47r 

7V4 ^2 



/(") - ^ 



1 



7 



187r 247r 647r 



In TV 



37r 27rn^ 

427 - 77 
3847r 



dn = — 
8 



71 







1 



37r 27rn2 



-0.02955870986828890 , 



lim 

N^oc 



1 

24 



1 

4! Jo 



n f(n)dn 



dn 

6 



Ar3 



/(") - 



1 



21 



37r 27rn2 32?™* 



dn - 



7Af 



907r 1447r 2567r 

dn = 



'dn 27rn2 32T:n^ 



128 



427 - 77 
3847r 



(A8) 



Note in particular how the first terms in the expansion of The complete series expansion of A is therefore obtained 
fin) act as counterterms, making all the integrals con- as 
verge, and the numerical calculation therefore straight- 
forward using Mathematica. The coefficient in front of 
the logarithmic terms is easily found from IjASp and ljA6p : 

4 1 1 

A = -t- 0.089237810/ir - —fir\n fir 

6nfir 3 27r 

+ -fi^r^ - 0.029558710mV + —fi^r^ \nfir 
8 647r 



E 



c 



2k 



\ (2fc - 1) 



ifir) 



2k-l 



2-K 



+ (A9) 

647r 



- Y^A*''^' + Oifi\'') . 



(AlO) 
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2. Large distances, ^ 1 

In the limit of large distances, the series expansion of 
A is obtained a lot easier than in the previous section, 
since we can now integrate the expansion of f{n) in (|A2|I 
for small n directly from to oo. The result is therefore 

2 16 -12 In 2 ilnfir 29.4398446730 



(120 In 2 



4„4 



' 204) In 60(ln^r)2 



-0(l//r8). 
(All) 



The error that follows from using the expansion (jA2l) all 
the way to n = cx) is of exactly the same order as what we 
would get from the next term in the expansion, which is 
easily shown: If we denote the first terms in the expansion 
by h{n), such that f{n) — h{n) + 0{n7) for small n, the 



relative correction to the gravitational potential can be 
written A = f{n)e-f''^"'dn = h{n)e- '''''' dn + 6, 
where the error 6 is bounded by 



l'^l< 



f{n) - h{n) 



dn . 



(A12) 



It turns out that f{n) — h{n) — (In rt)'^+ higher order 
terms. A rough estimate E{n) of the difference f{n) — 
h{n) which satisfies E{n) > \f{n) ~ h{n)\ for all n can 



therefore be chosen as E{n) = [l — (Inn)'^ 
means that 



This 



1^1 < 



'E{n)dn ■ 



(In 



higher order terms. 



(A13) 

i.e. precisely the same leading order term as what is 
obtained from the first term not included in IIA2II. 
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